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lordinate range and angle) so that spatial correlations and 
ipatterns can be identified. The format of these displays 
jis identical to that used for the simulation output. These 
idata (either real or simulated) are then subjected to 
ivarious processing procedures designed to enhance the 
Ipatterns and determine vortex locations. Each of the real 
lor simulated scans can be processed by one of three different 
procedures: simple frequency or wavenumber filtering, 

matched filtering, and deconvolution filtering. The final 
output is displayed as contour plots in an x-y coordinate 
I system, as well as in the form of vortex tracks deduced from 
the maxima of the processed data. 

In order that two LDV systems may be compared, data from 
either system is first converted to a common Cartesian 
coordinate system. This is done for each of the spectral 
moments in turn. Vortex tracks can either be obtained by 
cross correlation of data from both systems or by a triangu- 
lation method. A detailed analysis of run n\imber 1023 and 
run number 2023 is presented to demonstrate the data analysis 
procedure. Vortex tracks and system range resolutions are 
compared with theoretical predictions . 
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I. 


SUMMARY 


A data analysis program has been constructed to assess 
the SLDV (Scanning Laser Doppler Velocimeter) system per- 
formance, to validate the simulation model, and to test vari- 
ous vortex location algorithms. The program takes either real 
or simulated doppler spectra versus range and elevation, 
calculates the spatial distributions of various spectral 
moments (mean velocity, variance, skewness, kurtosis) or 

other spectral characteristics V I , , I 

^ peak' max' width' peak' sura 

and presents these in a 3-D display (coordinate range and 
angle) so that spatial correlations and patterns can be 
identified. The format of these displays is identical to 
that used for the simulation output. These data (either real 
or simulated) are then subjected to various processing pro- 
cedures designed to enhance the patterns and determine vor- 
tex locations. 

Each of the real or simulated scans can be processed 
by one of three different procedures: simple frequency or 

wavenumber filtering, matched filtering, and deconvolution 
filtering. The final output is displayed as contour plots 
in an x-y coordinate system, as well as in the form of vortex 
tracks deduced from the maxima of the processed data. 

In order that two LDV systems may be compared, data 
from either system is first converted to a common Cartesian 


coordinate system. This is done for each of the spectral 
moments in turn. Vortex tracks can either be obtained by 
cross correlation of data from both systems or by a triangu 
lation method. 

A detailed analysis of run number 1023 and run 
number 2023 is presented to demonstrate the data analysis 
procedure. Vortex tracks and system range resolutions are 
compared with theoretical predictions. 



II. INTRODUCTION 


The objectives of the data analysis study are: 

1) to identify the parameters affecting the detection 
capabilities of the LDV system, and 

2) to define and optimize vortex location algorithms. 
Analyses of a particular flight test (Runs 1023 and 2023, 
flyby of a B737) have been carried out. A detailed presen- 
tation of the results is given in this report. The format 
has been chosen to match that of the simulation code to 
facilitate detailed comparisons. Since vortex location is 

a prime objective, emphasis is placed upon the spectral 
distribution of the return signal and the effects of thres- 
holding (both of intensity and velocity) upon it. Vortex 
location algorithms using data from either one or two LDV 
systems are tested against the simulation results. 

The LDV return signal is in the form of a spectrum 
(intensity versus velocity) at a large number pf points in 
space. The positions of these points are given in terms of 
range R and elevation angle 9 for each system. Values of 
the following derived variables are computed for each spec- 
trum: 

1) average velocity, V 

2) variance, a 

3) skewness , 3 

' 3 ' 


4) kurtosis, k 

5) velocity (above threshold) at which the inten- 
sity (above threshold) is maximum, V 

6) maximum recorded velocity, 

7) difference between maximum and minimum non-zero 
velocity after thresholding, 

8) sum of the intensity in all velocity channels, 

^sum 

9) maximum intensity, 

These variables are plotted both as functions of range at 
constant elevation angle and as contour plots in Cartesian 
coordinates. The R-0 plot shows range-angle correlations, 
while the contour plot shows the spatial variation in the 
scan plane. 

The effects of applying intensity and velocity 
thresholds to the spectra have been examined. The optimal 
intensity threshold is found by successive approximation, 

increasing until a clear contrast is achieved between 
the localized vortex return and the distributed background 
return. The results of this process can be summarized as 
follows: 

1) A low intensity threshold results in a noisy "vortex" 
signature broadly distributed in space for all of the 
variables derived from the spectrum, except higher 
order moments of velocity which show localized signatures 
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2) Medium intensity thresholds result in I (and also 

■* sum 

^peak' non-zero velocity thresholds) being localized 

near the vortex core. Other variables are noisy and 
broadly distributed. 

3) High intensity thresholds result in localized distri- 
butions for all variables; I is most localized. 

sum 

Only two velocity thresholds (3.815 and 20.165 m/sec) 

have been examined, an.d no systematic investigations of the 
effects of the value of have been carried out. All 
variables (and especially are affected by 

The sensitivity of the results to the values of the 
thresholds is reduced by averaging the spectra. The averag- 
ing process redistributes the signal intensity in a given 
channel to its nearest neighbors according to a prescribed 
distribution formula. Averaging was found to be helpful in 
removing noise from either atmospheric turbulence or the 
LDV system itself. Three filtering procedures (simple 
frequency, matched filter, and deconvolution filtering) can 
also be used to enhance the maxima of the contour plots and 
to reduce the effects of noise. 


The simulation studies show that the core location 
is that of the maxima of V, a, 8, and k. Applying this 
criterion to the contour plots from the data tape gives 
the vortex core location for the flyby case. The vortex 
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track obtained by plotting the positions of the maxima in 

I ' derived from the data tape versus time can be analyzed 
sum 

to yield values of descent rate, vortex separation, and 
circulation. These values agree well with predicted values. 

Simultaneously obtained data from two LDV systems 
can also be used to determine the vortex track. Two dif- 


ferent procedures have been examined. In the first, the 
vortex location is prescribed as that of the maximum in the 
product of a derived variable for one system with that of 
the same variable for the second system (e.g. , 
the second, the vortex system is located by triangulation 
using the values of from both systems. The way in 

which information is used depends upon range. At 

long range, the elevation angle of the vortex core was shown 
by the simulation studies to coincide with that of the maxi- 
mum in At short range 70 m) , is a minimum 

at the vortex core, and maximum to either side in elevation 
angle. These criteria define a vortex core elevation angle 
for each system, and the core position is at the intersection 
of lines drawn from each system at the determined angle. 


6 



III. DATA REDUCTION AND ANALYSIS 


• Smoothing Procedure for the Signal Intensity Spectriim : 

A typical spectrijm is shown below: the abscissa is 

the velocity or channel numbers (the channel width is equal 
to ,545 m/sec) and the ordinate is signal intensity. The 
high velocity channels have little signals, the low velocity 
channels contain wind information, and all channels contain 
fluctuations. Even after thresholds are applied, there remain 
fluctuations which make it difficult to identify the vortex 
location. 

We have implemented a smoothing scheme which redis- 
tributes the intensity level I (n^> of the original data to 
its neighboring channels according to the relation 



The signals at channels below zero velocity are reflected 
back to channels above zero velocity. The result is a much 
smoother spectrum and contours obtained are also smoother. 
The horizontal line represents the intensity threshold 
The left figure shows the original spectrum, while the right 
figure shows the same spectrum after smoothing. 
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The variation is much smoother and gradual than the 
original spectrum. All variables derived from this 


smoothed spectrum should also exhibit similar smoothness. 



• Determination of Intensity Thresholds 

Since both the LDV system mechanisms and the atmos- 
pheric turbulence generate a certain amount of noise, the 
output from raw data without applications of intensity 
thresholds appear irregular and bear no similarity to the 
simulation results. It is therefore important to determine 
the noise level from the data. The noise level in general 
is a function of both channel number and time t. In princi- 
ple, one could block the laser beam window and obtain the 
j^ystem noise spectrum as a function of channel number at a , 
given time and use this in reducing the data for a subse- 
quent scan. In practice, this is not done; a constant 
noise level is specified for all channels. Applying this 
noise level to the data is straightforward. However, if 
the noise level is set too high, valuable signals will be 
lost, and if too low, the desired signature can still be 
indistinguishable from the background noise. A second 
method we have adopted is to look through the entire scan 
frame and select tho lowest signal intensity for each channel 
as the noise level. 

Different intensity threshold levels can be then 
applied above the noise level determined by either method to 
study the effects. 


9 



• Construction of the Contour Plots : 

Since data are only available at the focal points 
depicted by the fan beam configuration, it is convenient 
to redistribute the data on a rectangular mesh to obtain 
the contour plots. The separation between data points is 
large near the long range end of the fan, and it is desir- 
able to find a way which can fill up the void space between 
each range scan smoothly. 

The procedure we have implemented here is to dis- 
tribute the quantity carried by the focal point P (see 
sketch) onto a square having the same dimensions as that 
of the desired mesh. Then a weighted value (proportional 
to the shaded overlapping area) of that quantity is assigned 
to each of the neighboring cells A, B, C and D. The mesh 
dimension should be comparable or smaller than the largest 
separation between points in the fan beam. The void space 
outside the fan was assigned a value equal to the minimum 
of the quantity to be contour-plotted. This was found effec 
tive in removing the large number of contours which would 
otherwise appear near the edge of the fan. 




• Filtering Procedures : 

Techniques commonly used in the analysis of radar 
data to maximize the signal-to-noise ratio can also be applied 
to the LDV returns. Detailed descriptions are given by 
Skolnik (1970) . 

The return signal s (x) can be considered as consist- 
ing of two parts, f (x) which is due to the vortex at the 
unknown location x^ (which is to be determined) and n(x) 
which is due to noise. Then 

s(x) = cf (x-x^) + n(x) 


where c is an arbitrary constant. The Fourier transform 
of this equation is 



1K*X. 


cfg e 



Dividing both sides by fjv, we have 




ik*x 


ce 



In general, filtering involves applying a weighting function 
Wj^ to the Fourier components p^ and then retransforming to 
the physical plane. 
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F"^ (PjJWjj) 


= oj Wj 


xk»x -ik*x 
e ° e die + f’ 


1 /■'rev 

j 


The object of the filtering procedure is to enhance the first 
term with respect to the second. Different forms of the 
weighting or filter function Wj> produce transformed signals 
with different properties. 

Deconvolution filter ; Here we set Wjv = 1 # so that 




< 

I 


iicMx -x) ^ 

e ° dk + F <"k/*)c> 


= c6(x-x^) + G (x) 


In the absence of noise, this filter produces a delta 
function at x^. When substantial noise is present, G(x) 
can be appreciable at some values of x, and it becomes 
difficult to identify the delta function. 

Matched filter ; This is the optimum filter when 
the signal is observed in the presence of white noise 
(ng independent of k) . The weight function is wj^ - l^k^^ ' 
and 
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The underlined integral is a correlation and can be shown 

to have a maximum at x . Then the inverse transform 

o 

-1 

F will also be a maximum at x^. 

Simple frequency filter ; In this case, the weight 

y 2 

function is w^ = exp(-|k| ). This filter only smooths out 
the noise. 

The filters described above have been incorporated 
in the LDVDA code. The usefulness of filtering procedures 
has not been explored to any great extent and deserves 
further study . 

• Data Reduction Procedures ; 

The basic set-up is described in the Appendix, The 
data tapes from NASA/MSFC are first converted into 60-bit 
binary word tapes to be analyzed on the CDC 7600 computer. 
The data reduction program LDVDA is then applied to study 
these tapes. 

13 



In general, each range scan is composed of approximately 
50 points, and about 25 lines of sight for a complete angle 
sdan so that there are approximately 1250 range points in 
each frame. Assuming each focal point is separated by .008 
seconds in time, one complete scan will cover about 5 seconds. 
Each tape covers about 50 seconds of data so that there 
should be 12,500 data points altogether. For each focal 
point in space, there are six numbers to record its frame 
sequential number, x,y coordinates, range r, angle 6 and 
tl^e clock time plus one number to record the signal intensity 
for each of the 104 velocity channels. The computer cannot 
store all this information in the memory core at once, and 
therefore only 1295 data points are read and stored in the 
large memory core (LCM) at any one time* Since the record- 
ing of the data does not always start at the range maximum 
or minimum, portions of the 1295 data points are discarded 
before the beginning of the range scan is found, if the 
total number of data points in one complete scan frame to 
be read exceeds what is left in the LCM, the subroutine 
PUSHUP will discard data points which have already been used 

i 

aiid will read another 1295 points into the LCM. The noise 
level is found, and various intensity and velocity thresholds 
are applied to calculate V, o, g, k . 

^sum' ^peak- 

14 



IV. RESULTS AND DISCUSSION 

The data reduction program is capable of analyzing 
both flight and simulation data tapes. Thfe first output 
of both data tapes is the fan beam configurations. These 
are shown in Figures la and lb, where Figure la is the 
Run 1023 data and Figure lb is the simulation data. Both 
horizontal and vertical coordinates are of the same scale 
in meters. The caption near the top of the fan describes 
the lowest and the highest elevation angles, as well as 
the initial and final clock times in seconds. The jittering 
of the angle and range scan of the LDV system can be assessed 
from this plot. The simulation data (Figure lb) uses a fan 
of lines-of-sight separated by a constant angle of 1.5® and 
is extended from the aperture all the way to the outer • 

boundaries of the computation mesh. The higher density of 
range points near the vortex locations is a result of finer 
integration steps to resolve the higher velocity gradient 
there. 

Due to the large volume of output for each data tape 
and the large number of data tapes processed, we will con- 
centrate on discussion of only the best representative cases, 
i.e.. Runs 1023 and 2023, throughout this report. 

Run 1023 (run 2023 represents output from the second 
LDV system) was conducted at MSFC/NASA on July 20, 1974, at 

' ■ 15 . 






the Redstone Arsenal airstrip. A B737 was flown from 
north to south. The first LDV system was located at 61.11 
meters to the west, and the second system was at 240 meters 
to the east. 

Figure 2 shows four typical spectra . Figure 2a is 
a spectrum of the ambient wind; Figure 2b is that of the 
vortex. Notice the peak intensity is shifted toward higher 
velocity channels. When there is no noise, defining the 
noise as the minimum intensity in each channel should yield 
"noise” which is identically zero. Since there is no noise 
generated in the simulation data, this procedure did produce 
a zero noise level. Figure 2c shows a real spectrum due 
to ambient wind, and Figure 2d shows a spectrum due to the 
vortex, which brings nearly to 7 m/sec from the ambient 

wind of 1.63 m/sec. 

Figure 3 shows the applications of the noise threshold 
schemes to the real data spectra. The abscissa is the velocity 
in m/sec or is equivalent to the filter channel numbers. 

The ordinate is the intensity level. Figures 3a and 3b show 
the first intensity threshold scheme being applied to two 
typical vortex signatures, and Figure 3c and 3d show the second 
noise threshold scheme as applied to the same spectra. The 
dotted lines are the recorded spectra; the undotted lines are 
the noise. From Figure 3b, it is expected that the constant, 
intensity threshold of 51 is too low due to the spurious spike 

17 




















at a velocity equal to 18.64 m/sec. From Figures 3c and 3d, 
it can be seen that the second scheme should subtract out 
the ambient wind signature and leave the vortex signature 
intact. One might conclude that the second scheme should 
be a better scheme to achieve maximum clarity for the vortex . 
signature. However, as it will be shown later, if the 
threshold of the first scheme is high enough, the results 
are very close to that of the second scheme at lower inten- 
sity thresholds. 

Figure 4 shows the results of applying two different 

constant intensity thresholds to a set of spectra. Figure 4a 

shows a typical spectrum and a threshold level of 51. The 

values of V and V which result are shown for all R and 

wxdth 

0 in Figure 4b and 4c (left col\imn) . In the right column. 
Figures 4d, 4e, and 4f, the threshold is set at 60. The 
abscissa of each plot is range along the line-of-sight in 


meters. The ordinates are indicated by the captions; all 
have, units of m/sec. Each curve represents the variable 
versus the range, and each is displaced vertically by a 
constant distance proportional to the difference in elevation 
angle. The location of the LDV system, the lowest and high- 
est angle and the clock time recorded for the beginning and 
end of the frame are written on the plot. The minimum and 
maximum values of the variable appearing in the frame are also 
indicated. Vertical spacing between each curve represents 
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Figure 4, Effects of Different Constant Intensity Thresholds on V and V 
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half of the difference between the maximum and minimum of 

the variable plotted. Due to the appearance of several 

spikes at higher velocity channels in Figure 4a, it is 

expected that the V and in the left column will be 

less smooth than those in the right column. This is found 

when one compares Figure 4b to 4e and Figure 4c to 4f. 

Figures 4e and 4f show the localized bumps at elevation 

angles of .43 and .32 radian and at ranges near 109 meters 

and 127 meters, respectively. 

Figure 5 shows the effects of the same pair of 

thresholds on V and variance. The V shown 

peak' max' peak 

in Figure 5a is very noisy, while Figure 5d shows well 

localized vortex signatures. This indicates that is 

Very sensitive to the intensity thresholds. On the other 

hand, (Figures 5b and 5e) shows no dependence upon the 

intensity thresholds, as is expected. Therefore, we believe 

V should be a better quantity to work with than V , 

max j peak 

as far as the sensitivity to the intensity thresholds is con- 
cerned. Figures 5c and 5f show the effects of different 

2 

thresholds upon the variance a in (m/sec) . They show the 
same trend as appeared in but the degree of sensitivity 

is lessened. Figure 5f shows well localized patterns at 
angles of .43 and .31 radians and at ranges near 136 meters. 
From the simulation outputs which were discussed by Thomson 
and Meng(1974), it is known that if the line-of -sight passes 
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Figure 5. Effects of Different Constant Intensity Thresholds on 







directly through the vortex core there is no signal. One 

may ask whether the signatures in Figure 5f are due to one 

or two vortices. The separation between bumps at angles 

of .43 and .31 can be calculated to be approximately 13.6 

meters and the horizontal separation distance is 12.7 meters, 

which is shorter than the initial vortex separation ('^>18 

meters) predicted for the B737, but longer than the predicted 

vortex core diameter of four meters . Figures 6a‘ and 6b show 

the skewness (in (m/sec)^) and kurtosis (in (m/sec) at an 

intensity threshold of 51. They appear quite noisy. Figures 

6d and 6e show the skewness and kurtosis for = 60. Both 

th 

show well-localized maxima at an angle of .31 radians and 

have sharper angle resolution than that of the variance in 

Figure 5f. Figures 6c and 6f show almost identical results 

of the which is virtually independent of the value of 

the intensity threshold. has maxima near the vortices 

and varies smoothly with both angle and range. 

Figures 7a and 7b show This is not a good 

indicator for vortex location when there is no velocity 

threshold but, as will be shown later, it is a good indicator 

if a velocity threshold is applied. Figures 7c and 7d show 

i, . . , - . : 

the maximum in m/sec for each line-of-sight (that is, 

the highest for one complete range scan) versus the 

elevation angle in radians. Because of the low intensity 
threshold, the maximum in Figure 7c is as high as 

' . ' ■ ' 24'' 
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Figure 7. Effects of Different Constant Intensity Thresholds on I 
IIaximum vs. Angle 
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21 m/sec and the curve does not reveal the vortex location. 
Figure 7d shows two well localized peaks at angles of ,32 
and .43 radians. These angles are consistent with all 
previous plots. 

No velocity threshold is applied on the spectrum 
shown in Figure 8a. The vertical line indicates the zero 
velocity channel, the lower undotted line represents the 
^min determined by the second scheme, and the higher 

undotted line represents an intensity threshold of three 
added onto the level. Figure 8b shows the same spectrum 

with velocity threshold at 7.63 m/sec. The 'a' series of 
Figure 9 through Figure 17 are the results for "no velocity 
threshold" cases, and the 'b' series are those for the 
"velocity threshold" cases. 


Figure 9a shows V is broadly distributed. It becomes 

a minimum at vortex locations, while Figure 9b, in Which V 

is' more localized, shows the same trend. Figure 10a shows 

that V is well localized without a velocity' threshold, 
max 

and Figure 10b indicates that at this low intensity threshold 


level, velocity threshold produces a very noisy plot of V 


max 


Figures 11a and 11b are identical, since at low intensity 
thresholds will not depend upon the velocity thres- 
hold. Figures 12a and 12b is almost identical, 

and both figures show the same trend observed in V 
(Figure 9). Figures l3a and 13b show slightly improved 
localization for /F over that shown in Figure 9. The same 
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Figure 8, Spectra Shovmng Low Intensity Thresholds and 
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Figure 9, Effects of Velocity Threshold on V at Low Intensity 
Threshold 
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Figure 10. Effect of Velocity Threshold on at Low 
Intensity Threshold 
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Figure 11. Effect of Velocity Threshold on at Low 
Intensity Threshold 
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Figure 12. Effect of Velocity Threshold on at Lc 

Intensity Threshold 


Figure 13 . 


Effect of 
Intensit\ 



KANCr *10N6 TnI lINI Or SISNT IN N 

B 

Velocity Threshold on V variance at Low 
Threshold 
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Figure 15 , Effect of Velocity Threshold on Coefficient of 
Kurtosis at Low Intensity Threshold 
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Figure 16 . Effect of Velocity Threshold on at Low 
Intensity Threshold 
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Figure 17. Effect of Velocity Threshold on at Low 
Intensity Threshold 37 









trend that the velocity threshold brought about in Figure 9 
is again observed in Figure 13. 

The coefficient of skewness now as 

B = tI(v-7)^ l(v)lVa^ 

isj shown in Figure 14, This quantity is always positive. 
Its location of maxima coincides with the minima of a. 

Both Figures 14a and 14b show good vortex signatures, and 
velocity threshold further improves the localization of 
vortex signatures. The coefficient of kurtosis defined as 

K = [J(V“V)^ . I(v) 1/a^ 


is shown in Figure 15, Both Figures 15a and 15b show good 
signatures. The plots of skewness and kurtosis as defined 
previously appear quite different from those shown in 
Figures 14 and 15. Both skewness and kurtosis plots were 
fbund to be similar to the variance shown in Figure 13, but 
were quite noisy. The redefined coefficients of skewness 
and kurtosis seem to yield clearer signatures. 

Figure 16 shows that I _ is maximum at the vortex 
^ sum 

location, while the velocity threshold case in Figure 16b 
shows the opposite trend. Finally, Ipeak 
Figure 17. The velocity threshold is seen to be able to 
improve the angular resolution in significantly. The 

V , versus angle plot for this low intensity threshold case 

p03.JC 


is not informative since is equal to the constant velo- 

city at the highest channel for all angles. 
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Figure 18 shows the effects of different variable 

I 

intensity thresholds above the noise level, determined by 
the second scheme. Figures 18a and 18d show the recorded 
spectrum as the dotted line. The undotted line is the 
threshold spectrum obtained by adding a constant value 
(20 or 25) to the noise level derived from the minimum in- 
tensity in each level. The threshold spectrum shown was used 
in analyzing the data; although the data spectrum shown lies 
mostly below the threshold spectrxam, other spectra do not, 
and these are responsible for the signals. One might expect 
to find the same effect as was found in the first scheme. 


However, the average velocity in Figure 18b shows a single 
peak at a location totally unrelated to the vortex locations. 
This was because of a very high intensity return at that 
range point, which was obviously a noise return. Due to 
that very high signal, all the rest of the variation is 
suppressed. Figure 18e shows maxima in V, but it is hard 
to identify a localized vortex signature. Figures 18c and 


18f, showing 3^^^© almost identical, 

noise point which appeared in the V plot. 


except for the 
Because of the 


lower intensity threshold. Figure 19a shows several peaks 


appearing in the plot of while Figure 19d shows only 

one peak. Figures 19b and 19e show and Figure 19e 

appears to be similar to Figure 5e. Figures 19c and 19f, 


together with 18b and 18e, show that plots of the variance 





Figure 18 , Effector Different Intensity Thresholds Above Doise Level 
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a and of V depend in a similar fashion upon the intensity 
threshold. Figures 20a and 20c, and 20b and 20d, show 
that skewness and kurtosis also behave similarly. Figures 


20e and 20f show the maximum versus elevation angle. 

Figure 20e has the higher maximum value of and shows 

three peaks at angles of .25, .35, and .45 radians. (These 
are also found in Figure 19a.) Figure 20f shows the highest 


''peak 

and .40 radians. 


with three peaks at angles of 




Results from applying a higher intensity threshold 
are shown in the a, b, and c portions of Figures 21 through 
24 . These plots exhibit a behavior very similar to the 
series shown in the d, e, and f portions of Figures 4 thrbugh 
7. The second noise determination scheme (minimum intensity) 
does not produce clearer vortex signatures than the first 
scheme (constant noise level) , and requires more effort to 
apply. However, the second scheme should be less sensitive 
to arbitrary assumptions in setting the intensity threshold* 
Figures 21 through 24 show the effects of velocity 
thresholds. The vertical lines (at V = 0 in Figure 21a 
and V = 7.3 m/sec in Figure 21d) show the threshold values. 

A constant value of 30 is added to the minimum-signal noise 
level in all cases. All signals below the velocity threshold 
are dropped. This will subtract all wind information for 
the right column plots, so that signals will be due to the 
vortex only. 
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Figures 21e and 21f show signals which are better 
isolated than those in Figures 21b and 21c. Figures 22d 
and 22e show almost identical results as those in Figures 
21e and 21f. However , Figures 22f and 23 d,e,f reveal 
only one peak instead of the four peaks which appeared 
in the (Figure 22d) and (Figure 21f) plots. 

The most remarkable effect of a velocity threshold is on 
^sum 23c and f ) and (Figures 24a and b) ; 

Ipeak ^ good indicator of vortex location if velocity 
thresholding is applied. Figures 24c and 24d both show 
two peaks at angles of .30 and .43, and the velocity 
magnitudes are the same. The gradual drop off on both 
sides of the peak shown in Figure 24c disappears when 
velocity thresholding is applied. 

Figures 25 through 32 show a series of contour plots 

W-'^peak' ^sum' ^peak’ ^ constant inten- 

sity threshold of 60 is applied and there is no velocity 
threshold. 


The method of preparing the contour plots was dis- 
cussed in Section III, The grid size is 200 m x 200 m. 
Values of the quantity on the contours are indicated near 
the top of each plot. An auxiliary cross-sectional plot 
is given below each contour plot to demonstrate the vari- 


ation versus both x and y distances. The time of the scan 


is recorded on both the contour and the cross-sectional 
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Figure 26 . Contour Plot for 
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Figure 27. Contour Plot for V 
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Figure 28. Contour Plot for a 
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Figure 29 . Contour Plot for e 
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Figure 32 . Contour Plot for 
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plots. The quantity being plotted is indicated near the 
upper left corner of each contour plot. The highest and 
lowest value of the quantity are written out in the cross- 
sectional plots . 

Figure 25 shows contours of V in meters/sec. The 
separation distance between the two maxima in the contour 
plot is approximately 25 meters, which is quite close to 
the theoretical prediction of 28 meters. The vertical 
dimension of each maximum is about 25 meters, which is 
also reasonable when compared to the simulation results. 

Figure 26 shows the contour plot for V , It is quite 
similar to that of V, and the maxima are closely correlated 
to those of V. Figure 27 shows the contour plot for 
Due to the high intensity threshold, it seems that the 
left vortex signature has been lost in this plot. All higher 
moments of the spectrum seem to give well localized signatures, 
but they also lose the weaker signature at high intensity 
thresholds. 

2 

Figure 28 shows the variance a in (m/sec) . It shows 
only one peak at (125, 40); the other peak is almost com- 
pletely lost. Similarly, in Figxires 29 and 30, the contours 
for 3 and k are shown to have one peak only. Finally, 

Figures 31 and 32 show the contours of 1^^ and Ipgak’ ^sum 
shows well-defined maxima which correlate with those in V 

and V but disagree with the position of the maxima of 
max 
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all other quantities. The reason for the discrepancy is 
unknown and deserves further study. Locations of the 
maxima in (Figure 32) are difficult to determine 

unambiguously; they become fewer and more sharply defined 
when a velocity threshold is applied. 

Summarizing the discussion so far, we can conclude 
that the results from both noise level determination schemes 
appear to be similar. The choice of the intensity thresholds 
is very important for all variables, but it is less critical 

^pealc- ® locator, 

is a good vortex locator if a velocity threshold is applied. 
At high intensity thresholds, the low order velocity moments 
yield spatially-distributed, but smooth, contours. The 
high order velocity moments are very localized but lose 
the weaker signals completely. 


A series of contours for I „ are shown in Figure 

33abcd. It is known from simulation studies that maxima of 

V, a, g, and k are identified with the vortex location if 

range is above 100 meters. On the other hand, the data 

output shows the locations of maxima in V are closely 

correlated to those for 1^^^. Therefore, one could identify 

a maximum in I as the vortex core. Plotting the maxima 
sum 

of I as functions of time, we can obtain the vortex tracks 
sum 

depicted in Figure 34. It can be seen that there are two 
vortices moving from left to right at a rate of about 1.84 m/sec 

5g 












and descending at 'V'1.08 m/sec. The average vortex separa- 
tion is about 30 meters, so that the total circulation 

. 2 2 
contained in the vortex is about 204 m /sec (2205 ft /sec) , 

which is also close to the theoretical prediction of 2360 

ft /sec. The entire content of a seven track data tape 

can be represented by these nine points . 

The actual range resolution of the LDV system can 
also be assessed from the contour plots of I in Figure 
33. If we define the range resolution As as the half width 
of along the beam, we can determine the dependence of 

As upon tange s. Figure 35 shows As (meters) versus range 
(meters) obtained this way. Also shown is the predicted 
value As = 2f ^ where f = 2.1, X = 10. 6y and R = 15.24 cm. 

The range resolution is poorer than predicted value. 

Vortex tracks can also be obtained by analyzing 
data from two LDV systems scanning simultaneously over 
generally the same area. One procedure is to multiply the 
variable defining the vortex track (say V) from one system 
by its value given by the other system, and to then plot 
contours of this product. The positions of the maxima are 
then identified as the vortex positions. A second procedure 
is to apply triangulation, as discussed in Section II. 

Figure 36 shows an overlapped fan beam configuration. 

The first system (left) is located at -61.11 meters; the 
second system (right) is at 240.33 meters. The contour 
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Range Resolution in Meters 



Figure 35 . Range Resolution as a Function of Range 
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Figure 36. Fan Beam Configurations from Tv^fo LDV Systems 
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plots for V from both systems are plotted separately in 
Figure 37a and b. Both show elongation characteristic 
of the LDV system. The contour for their product is given 
in Figure 37c, which is well defined over a reasonable 
dimension, and the noisy contours elongated along the range 
appearing in Figure 37a and b are eliminated. Doing this 
for each time step will produce the desired vortex track. 

If the same process as described above is repeated 
for the quantity I _ , smoother and clearer vortex tracks 
can be expected. This was not done, but simple comparisons 
• ^sum outputs of these LDV systems were performed. 

The contour plots of I from both Runs 1023 and 

2023 at approximately the same time are shown in Figures 
38 and 39. The intensity threshold was set at 50 in 
Figures 38a and 39a, and at 25 in Figures 38b and 39b. The 
run numbers are indicated in each plot; ground level is 
marked in Figures 38a and 39a. The cross symbols are used 
to identify the maxima in each plot. Range in meters is 
measured from each LDV system and can be converted to a 
cpmmon frame. Figure 38a shows two maxima, while Figure 38b 
shows only one. The location of maxima when converted to 
the common frame are not closely correlated. Figure 39a 
shows the maxima at (53,35) and (76,25), while Figure 39 b 
shows the maxima at (70,25) and (88,27) when converted to 
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a common frame. There are considerable deviations. The 
discrepancy could be attributed to range resolution and 
different thresholding conditions. 

A series of plots of vs. angle for both LDV systems 

are given in Figures 40 through 42. The 'a' series is from 
the first system; the 'b' series is from the second system. 

From each of these plots, the angle at which the is 

a maximum can be read. Then a ray can be drawn at that 
angle from the location of the corresponding system. Inter- 
section of the rays will determine the vortex locations. 

The accuracy of the vortex track determined this way in 
Figure 43 is not very clehr. The only guidance is to com- 
pare the simulation results of versus angle to actual 

vortex location. It appears difficult to determine the 
vortex track from information alone. 
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Figure ^3. Vortex Track Determined by 
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V. CONCLUSIONS AND RECOMMENDATIONS 


Two vortex detection modes were considered: a 

"soft target" mode which deduces the vortex track by assuming 
the vortex is located at the maxima of a well-defined, smoothly 
varying quantity such as I_,,_; a "hard target" mode in which 
a triangulation method was applied to locate the center of 
the peak velocity for both LDV systems. Several conclusions 
and recommendations were drawn from examination of the data: 

1) The two noise-determining schemes are successful and 
both produce clear vortex signatures; the best value 
of I^j^ must be evaluated by trial and error. 

2) The vortex wake appears as a spatially-distributed 

target; the clarity of the signature depends critically 

on the chosen thresholds, either I^, or V^, . 

th th 

3) The spatial dimensions of V, V .appear to be 

larger than those for cj, 3, and k. However, the utility 
of these variables in locating the vortex core deserves 
further study. 

4) The variable which has the least dependence upon the 

value selected for I., is I and the maxima of its 

th sum 

contour appear to define the vortex center satisfactorily. 

5) Range resolution based on of the LDV system appears 

to be poorer than '“theoretical prediction. 
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6) The best versus angle plot is generated by the 

data analysis code when no velocity threshold is 
applied and only medium is applied. Too high an 

will eliminate the weaker vortex signal and too 
low an will yield unrealistically high If 

a velocity threshold is applied, the angle dependence 
near a vortex may be lost. 

7) Certain ambiguities must be cleared up before the 

triangulation method can produce reliable vortex tracks. 
This is mainly because the versus angle plot is 

range dependent; that is, it appears very different 

at long range from that at close range. In addition, 
it requires better coordination between the two LDV 
systems to maximize the overlap in scanning. 

8) The contour plot of the cross correlation of I.,,- from 

tv/o LDV systems should be a good vortex track generator. 

Evaluation of I should also be easily adaptable in 

sum 

the present on-line system to yield the vortex track. 

9) The noise level as a function of time should be recorded 
on the high speed tape to eliminate the trial-and-error 
procedure. 

10) A translator should be included in the LDV system to 
obtain the sense of velocity. The velocity vector 
obtained from the existing data is not meaningful. 
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Appendix; FLOW CHART FOR THE LDV DATA ANALYSIS PROGRAM 


The basic program setup is described in the flow 
chart. The data tapes are first converted into 60-bit 
binary words in order to be analyzed on the CDC 7600 
computer. The data reduction program LDVDA is then applied 
to study these tapes. 

Since the recording of the data does not always 
start at the range maximum or minimum, the first range 
maximum or minimum is used as the first point for data 
processing. This is done in the subroutine RESET. Then 
the stored data points are fetched and the spectrum plotted 
for several space points. All the variables v, 

Vak' W' ^sum’ ’^peak ^hen computed as a 

function of ranga and elevation angle. 

After one complete scan is obtained, the subroutine 
DISPLAY is called to present various plots. The spatial 
pattern of data points is first displayed. The subroutine 
3DPR0FL is called to produce the r,6 plots for all the 
variablest v, 0. B, k, 

The highest velocity for each angle is also plotted versus 
angle. 


The filtering procedures are then applied to analyze 
the data. First of all, the fan beam configuration is mapped 
into a rectangular mesh. This is done by calling the sub- 
routine SETUP; then CONMAP is called to generate the contour 
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plots for the array on the rectangular mesh. This is 
valuable for comparing to the plots after the filtering 
is! applied. Subroutine MATCH then performs one of three 
different filtering operations: matched filter, smoothing 

or deconvolution. . 

Subroutine CONMAP is again called to plot the 
contours of the filtered data and all the variables v, 

W' ''width' ^SUKI' Wk processed 

through the same contour plotting procedure , 



II. Flow Chart for the LDV Data Analysis Program 

The basic program setup is described in the flow 
chart. The data tapes are first converted into 60-bit 
binary words in order to be analyzed on the CDC 7600 
computer. The data reduction program LDVDA is then applied 
to study these tapes. 

Since the recording of the data does not always 
start at the range maximum or minimum, the first range 
maximum or minimum is used as the first point for data 
processing. This is done in the subroutine RESET. Then 
the stored data points are fetched and the spectrum plotted 
for several space points. All the variables v, 

Vak' %ax' ^sum' ^peak computed as a 

function of range and elevation angle. 

After one complete scan is obtained, the subroutine 

DISPLAY is called to present various plots. The spatial 

pattern of data points is first displayed. The subroutine 

3DPR0FL is called to produce the r,0 plots for all the 

variables: v, v v ,rV /a/B/K/I #1 

V, max - ' sum' peak 

The highest velocity for each angle is also plotted versus 
angle. 

The filtering procedures are then applied to analyze 
the data. First of all, the fan beam configuration is mapped 
into a rectangular mesh. This is done by calling the sub- 
routine SETUP? then CONMAP is called to generate the contour 



plots for the array on the rectangular mesh. This is 
valuable for comparing to the plots after the filtering 
is applied. Subroutine MATCH then performs one of three 
different filtering operations; matched filter, smoothing 
or deconvolution. 

Subroutine CONMAP is again called to plot the 
contours of the filtered data and all the variables v, 

max' wxdth' ' ' sum' peak ^ t' 

through the same contour plotting procedure. 




SETUP 

;to seed the fan beam 
jfocal points on a 
;flxed x,y-plane mesh 


I CORMP 

generate contour plot 

r~ 

MATCH 


PUSHUP " 

discard used data 
points and read addi- 
’tional data to fill up' 
1295“point LCM storage 



CALL CONMAP : 

1 

- ■ ■ ----- f. 

CALL MATCH \ 

i 

apply the matched filter ; 
techniques 

CALL CONMAP 
CALL PUSHUP 

push out used data and 
continue to read data 





















